Fission-fragment mass distributions from strongly damped shape evolution 
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Random walks on five-dimensional potential-energy surfaces were recently found to yield fission- 
fragment mass distributions that are in remarkable agreement with experimental data. Within 
the framework of the Smoluchowski equation of motion, which is appropriate for highly dissipative 
evolutions, we discuss the physical justification for that treatment and investigate the sensitivity of 
the resulting mass yields to a variety of model ingredients, including in particular the dimensionality 
and discretization of the shape space and the structure of the dissipation tensor. The mass yields are 
found to be relatively robust, suggesting that the simple random walk presents a useful calculational 
tool. Quantitatively refined results can be obtained by including physically plausible forms of the 
dissipation, which amounts to simulating the Brownian shape motion in an anisotropic medium. 

PACS numbers: 25.85.-w, 24.10.-i, 24.60.Ky, 24.10.Lx 



I. INTRODUCTION 

A recent study performed random walks on five- 
dimensional potential-energy surfaces and extracted 
fission-fragment mass distributions that are in remark- 
able agreement with experimental data (l|. The striking 
simplicity of the calculation, its unprecedented predictive 
power, and the availability of tabulated potential-energy 
surfaces for essentially all nuclei of potential interest [2] 
raise the prospect that the method may provide a quan- 
titatively useful calculational tool for obtaining approxi- 
mate fission mass yields for a large region of the nuclear 
chart. 

However, before such applications can be made with 
any confidence, a number of issues need to be clarified, 
in particular why such a simple, and somewhat arbitrary, 
treatment can yield such apparently good results. Be- 
cause of the potential utility of the treatment, we address 
these issues below. 

As discussed already in the pioneering papers by Meit- 
ner and Frisch [3j and Bohr and Wheeler 0, Q in 1939, 
nuclear fission can be viewed qualitatively as an evolu- 
tion of the nuclear shape from that of a single compound 
nucleus to two receding fragments. The character of the 
shape dynamics is still not well established but, as a step 
away from a purely statistical approach toward a full dy- 
namical treatment, it is interesting to explore scenarios 
in which the shape evolution is strongly dissipative. 

Early studies of fission dynamics with the so-called 
one-body dissipation suggested that the nuclear shape 
motion is strongly damped [6J, |7| and it was advocated 
that a reasonable starting point for determining the aver- 
age evolution would be to balance the conservative force 
provided by the potential energy with just the friction 
force resulting from the dissipative coupling between the 
deforming nuclear surface and the nucleon gas [3, Q ■ A 
simple feature of strongly damped motion is that the in- 
ertial forces are relatively unimportant so it is less crucial 
to know the inertial-mass tensor for the shape motion. 



The earliest numerical studies of dissipation in fission 
dynamics were focused on the damping effect on the mean 
motion only, using various physical models for dissipa- 
tion, including incrtial effects, and using macroscopic po- 
tentials [6j, 13, 13, LLQJ ■ Consideration of the stochastic force 
in fission dynamics began as early as 1940, when Kramers 
considered the average delay in establishing a stationary 
flow rate over a one-dimensional barrier, thus inferring an 
increase of the fission lifetime due to dissipation [llj . Fur- 
ther approximate treatments of the Fokker-Planck equa- 
tion in one or two dimensions began around 1980 [12| - 
[16J . These calculations often retained the assumption of 
a constant inertia and dissipation, with very simplified 
potentials. About a decade later, numerical investiga- 
tions of Langevin equations for dynamics including iner- 
tia, damping, and Markovian stochastic forces were be- 
gun by sever al g roups; reviews of such work were given 
in Refs. [2l|, [22[. These types of simulations continue 
[X7H19I ] . They employ two or three shape degrees of free- 
dom, macroscopic potential energies, fluid dynamic iner- 
tias, and more recent calculations usually use some form 
of one-body dissipation. Because of the use of macro- 
scopic energies, they are often explictly characterized as 
applying to systems with significant excitation energy. 

In contrast, the present investigation focuses on sys- 
tems having relatively low excitation energy, such as 
those produced by thermal neutrons, where it is essential 
to include microscopic (shell) effects in the potential en- 
ergy. For 5,254 even-even nuclei with 170 < A < 330, 
such potential energies have been calculated with the 
macroscopic-microscopic method on a five-dimensional 
lattice of 5,315,625 shapes Q. These potential-energy 
surfaces are the most comprehensive currently available 
and form the basis for our present studies. 



II. FORMAL FRAMEWORK 

We picture the nuclear fission process as an evolution 
of the nuclear shape from a relatively compact mononu- 
cleus to a dinuclear configuration. The nuclear shape is 
described by a set of parameters, x — {Xi}> whose time 
development is the result of a complicated interplay be- 
tween a variety of effects, as we now discuss. 

Most basic, and also most easy to understand, is the 
potential energy associated with a given shape, U(x), for 
which a number of relatively mature models have been 
developed. We employ potential-energy surfaces that 
have been calculated with the macroscopic-microscopic 
method in which the potential energy is the sum of shape- 
dependent macroscopic (liquid-drop type) terms and a 
microscopic correction that reflects the structure of the 
single-particle levels in the effective potential associated 
with the specified nuclear shape [201 ] ■ The potential en- 
ergy provides the driving force, F pot (x), which has the 
components F pot = —dU/dxi- 

The driving force seeks to change the nuclear shape 
and the associated matter rearrangement gives rise to a 
collective kinetic energy. Furthermore, the shape degrees 
of freedom are coupled dissipatively to the internal nu- 
clear degrees of freedom and, as a result, the shape evo- 
lution is both damped and diffusive. The nuclear shape 
dynamics should therefore be treated by transport meth- 
ods that allow for the stochastic elements of the dynam- 
ics. Stochastic transport approaches to nuclear dynamics 
have been reviewed in Refs. [2lll22|. 

The Lagrangian associated with the shape degrees of 
freedom has a standard form, 



c(x,x) = hJ2** M vbd*s- u bd 



and the (generalized) momentum can be obtained as 
tv = dC/dx, with the components being ni = Y] - MijXj- 
The inertial-mass tensor, M{x) 1 is not as well under- 
stood as the potential energy and cannot yet be calcu- 
lated with comparable accuracy. In the present investi- 
gation, the inertia is ignored because it is expected to 
play a relatively minor role if the dissipation is strong 
(see below). 

To a first approximation, the average effect of the cou- 
pling between the shape and the residual system is a sim- 
ple friction described by a Rayleigh dissipation function, 



Hx,x) = 



§Sxi7ti(x)*i = \Q i 



(2) 



which is equal to half the average rate at which energy 
is being transferred to the internal degrees of freedom. 
The dissipation tensor f(x) governs the associated fric- 
tion force F Mc ( X ), with f* ic = -dF/3 Xl = - £, laXj- 
We shall invoke the simple wall formula 7] to obtain es- 
timates of j(x) ( see Sect. IIII|> . 



The average shape evolution x(t) 1S then governed by 
the Lagrange-Rayleigh equation, 



dt 



dC 
dx 



OF 



(3) 



where the derivatives should be evaluated along the mean 
trajectory. When the damping is strong, the resulting 
motion is slow. In that scenario, the acceleration terms 
as well as terms of second order in the velocities may be 
neglected. The corresponding creeping evolution is then 
determined by the demand that the friction force exactly 
counter balance the driving force, i.e. F pot + F frlc = 0. 
This equation can then be solved for the average velocity 
X: Xi = Sj ^ijFf > where the mobility tensor /x is the 
inverse of the dissipation tensor 7. Because the dissipa- 
tion rate Q is always positive, the tensor 7 is positive 
definite; so its inverse /j, always exists and its eigenvalues 
are positive. 

The friction force represents the average of the inter- 
actions of the shape with the internal degrees of freedom. 
As is common, we shall assume that the remaining inter- 
action is stochastic and the associated force is denoted 
by F ran (t). (Its time dependence is indicated explicitly 
because it is expected to vary rapidly on the time scale of 
the shape evolution.) By definition, it vanishes on the av- 
erage, (F ran (t)) = 0, and we assume that its time depen- 
dence is Markovian, so (F/ an («)i^ an (t / )) = 2Ty bJ 5{t-t r ), 
where T(x) is the shape-dependent nuclear temperature 
(see later). 

The actual shape evolution is thus both damped and 
diffusive and the trajectory x{t) ls governed by the 
Smoluchowski equation of motion in which the driving 
force from the potential is counterbalanced by the full 



dissipative force, i.e. _F pot + F 



flic 



= 0. This con- 



W dition immediately yields the instantaneous velocity, 



X (t) = n-(F» ot + F iaa (t)) 



(4) 



We assume that the nuclear shape evolution is described 
by this equation and it is thus akin to Brownian motion. 
The net displacement accumulated in the course of a brief 
time interval At is then 



$X 



A/ 



X(t) dt = n 



F pot At 



A/ 



\t)dt 



(5) 



where we have chosen At to be so small that both the 
driving force and the mobility tensor can be considered 
as constant. The first term in the above expression ([5]) 
is deterministic and represents the average displacement, 
corresponding to the mean trajectory provided by the 
Lagrange-Rayleigh equation ([3]), whereas the second term 
is stochastic and arises from the inherent thermal fluctu- 
ations that give the evolution a diffusive character. 

Because of the diffusive nature of the dynamics, it is 
appropriate to describe the system by a probability den- 
sity P(x\ t). The average shape is then described by 



/ 



x(t) = / d X xP(x;t) 



(6) 



while the fluctuations around the average are character- 
ized by the correlation tensor cr(t) having the elements 

*•«(*) = / dx (xi - Ut)) (Xj - Xj(t)) P(x; t) . (7) 

When starting from a given shape specified by % , we 
have P(x>t = 0) = S(x — Xo) an( i the distribution then 
shifts and broadens in the course of time. In the Fokker- 
Planck approximation, this evolution is described by 



d_ 

dt J 



P 






P(x;t), (8) 



where V{x) is the drift coefficient (vector) and "D(x) is 
the diffusion coefficient (tensor). These transport coef- 
ficients are simply related to the early evolution of an 
initially narrowly defined distribution, namely 



|x = V = /x-FP°\ |<r = 22> = 2 M T 



(9) 



The drift rate x follows immediately from ([5]) and the 
diffusion rate <x is also readily obtained because the co- 
variance matrix for the changes in x is given by 

pAt pAt 

(SxSx) = / / W)x(t'))dtdt' (10) 

Jo Jo 

pAt pAt 

= /x-(F ran (i)F ran (t'))-At^di' (11) 

Jo Jo 

pAt pAt 

= / 2TfS(t-t')dtdt' = 27>At. (12) 

Jo Jo 

We note that VT = V ■ F pot in accordance with the 



fluctuation-dissipation theorem 23 1. 



A. Direct simulation 

One way to proceed is to solve the Fokker-Planck trans- 
port equation (JH)) which is a partial differential equation 
for a time-dependent function of N variables (in our case 
JV=5, see later). This is a formidable task and we instead 
perform direct Monte-Carlo simulations of the Smolu- 
chowski equation of motion (j4j to generate suitably large 
samples of individual stochastic shape evolutions. 

For this purpose, it is convenient to write the mobility 
tensor explicitly in terms of its eigenvectors {% }, nor- 
malized such that x ' X is the eigenvalue /i„ (which 
is always positive, as explained above), 



thj 



£x1 n) xf 



(13) 



Assuming now that the current shape is characterized by 
the value %, we wish to propagate the shape forward to 
a slightly later time t + At. The average shape change 
is readily obtained from Eq. §5§ as Sx = H • F pot At. 



The random contribution to the shape change is most 
easily sampled in the eigcnframe of the mobility tensor 
because the increments in each principal direction may 
be sampled separately. Invoking the eigen representation 
of (a (1131) . we may then obtain the total increment in x 
accumulated in the course of the small time interval At, 



Sx = E* W Aix (n) -* ,pot + V2TAl£r l 



(14) 



where {£n} are random numbers sampled from a stan- 
dard normal distribution having zero mean and unit vari- 
ance. This propagation procedure is easily implemented, 
once the mobility tensor fj, has been diagonalized. The 
average of the accumulated change is then 



8 Xl = AiVif'iW.F"', 






(15) 



because (£„) = 0, and the accumulated correlation cry- 
becomes 



a %3 = 2rA*£« B) *r. 



(16) 



because (£«£«.') = S nn >. Both are proportional to At. 

It is an important feature of the propagation scheme 
(fi"4f that when the time interval At is sufficiently small 
the generated ensemble of dynamical histories remains 
unaffected by a subdivision of the time interval, so the 
numerical solution of the transport problem is robust. To 
see this, imagine that the time interval used in Eq. (J14I) is 
subdivided into a number of shorter intervals, {At*-™'}, 
with ^2 m Ai^™) = At. If (i, F, and T remain unchanged 
during At, the resulting combined change in x becomes 

sx = E„^x (m) with 



m) 



8 x {m) = J2 X (n) &t {m) X {n) ■ F pot + v^TAiM & 

n 

(17) 
Recalling that (£„ ') vanishes, we see that the combined 
average change remains the same as before, 

$>X (m) ) -^x ( " ) (^A^) ) ^(«). F pot ; (18) 



as docs the accumulated covariance <x 



a- 



<(£*xr)(£*xr)> = 2r (E A* (m) )E 



~(n)~(n) 
v . V ■ 

At Ai 



because (£lT ^™ ') =S mm >5 nn 



(19) 

n £„™ ) =Smm>5nn>- Thus the diffusion pro- 
cess is robust under changes in the employed time inter- 
val At, as it should be. Of course, this invariance per- 
tains to the distribution of dynamical histories, P(x',t), 
whereas any individual trajectory does change when At 
is changed. 

A related invariance holds when the overall magnitude 
of the dissipation tensor 7 is changed: If the elements 
lijix) are a U increased by the common factor N(x) then 



the local time evolution proceeds at a rate that is N(x) 
times slower, but the resulting distribution of shape tra- 
jectories, P(x': t), remains the same. This convenient fea- 
ture allows us to arbitrarily rescale the friction locally 
to facilitate the numerical treatment, because we are not 
here interested in the actual time evolution but merely in 
the final distribution of mass divisions. Such local rescal- 
ing of the dissipation tensor is equivalent to adjusting the 
local clock rate which will obviously not affect the out- 
come of the process but merely changes how much time 
is spent at various locations. For convenience, we shall 
therefore assume that the eigenvalues Jnix) are one on 
average, for each particular shape X- 



B. Discrete random walk 

It is instructive to consider the simple situation when 
the mobility tensor is aligned with the lattice, i.e. fj, is di- 
agonal, fXij = Hidij. The transport process then reduces 
to a standard random walk, i.e. Eq. (TT4"f reduces to 



5 X > = luFiAt + y/2TiiiAt£i 



(20) 



If the potential energy U is known for any value of the 
shape parameter x the local force F can be obtained as 
the corresponding gradient, Ff° = —dU/dxi, and the 
transport process (f2T)]) can be readily simulated to yield 
an ensemble of evolutions. Because each increment 5\i 
is a real number (i.e. not necessarily an integer), each 
evolution is represented by a sequence of shapes whose 
coordinates {xi} may take on any fractional value within 
the considered parameter domain. 

However, the potential energy employed |2| is avail- 
able only on a discrete lattice whose 5,315,625 sites are 
labeled by the integers {IJKLM}, corresponding to inte- 
ger values of the shape coordinates, {xi, • ■ • , Xs}- By per- 
forming a pentalinear interpolation (see AppendixEJ, we 
may obtain an approximate representation of the poten- 
tial energy for arbitrary (fractional) values of the shape 
coordinates, U(x), an d the above continuous transport 
process (|20p may then be simulated. 

A simpler but approximate treatment of the above con- 
tinuous transport problem consists in performing a ran- 
dom walk on the discrete lattice, i.e. the shape parame- 
ters {Xi} take on only integer values. This can be conve- 
niently accomplished by means of the standard Metropo- 
lis sampling procedure (see below), as was originally done 
in Ref. |l| . The quality of this approximation depends on 
the lattice spacing and if the spacing is reduced then the 
discrete random walk becomes a better approximation to 
the continuous transport process. 

To understand how the continuous random walk (|20j) 
can be approximated by a Metropolis procedure, we first 
note that when the mobility tensor is diagonal, each lat- 
tice direction i is an eigen direction and can, therefore, 
be sampled independently. In the discrete treatment, the 
step size in the direction i is fixed by the lattice spacing 
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FIG. 1: (Color online.) The effect of the discretization of the 
shape space on the calculated fragment charge distribution is 
illustrated by comparing the results of the Metropolis walk 
introduced in Ref. J| with those of the corresponding contin- 
uous process; the experimental data [25|, [2(| are also shown. 
Throughout, as in Ref [1(, the calculated mass yields Y(At) 
have been transformed to charge yields by a simple scaling, 
Y(Z{) = (Ao / Zo)Y (A{) , where Zq and Aq are the charge and 
mass numbers of the fissioning nucleus. 



Aj, and we need to know the probabilities for taking a 
forward or backward step during a brief time interval At, 



>« 



V ( £At. 



The associated Fokker-Planck transport 



coefficients, which express the rate at which the mean 
location changes and half the rate at which the variance 
in the location grows, are therefore given by 



Vi 

V, 



(i/« - ^ } )Ai = ViFi , 



(21) 
I^f+^A^T, (22) 



where Fi is the force in the lattice direction i, and 
Hi = l/"fi is the mobility in that direction; they are in- 
dependent of the lattice spacing Aj. These relations can 
be readily solved for the rates, 



(i) _ Hi ,rp , 1 p a i Mi 



'TThAUi], (23) 



where AUi = — FiAi is the change in the potential as- 
sociated with an increase of x-i by A». It then follows 
that Px'/P: ~ ex p(~ AUi/T), which is precisely what 
characterizes the Metropolis procedure, as we now show. 
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FIG. 2: (Color online.) The effect of changing the lattice 
spacing for the Metropolis walk: The standard lattice Q] is 
modified in the Qi direction by introducing two additional 
lattice sites between the original lattice sites, thus tripling 
the density of lattice sites in that direction; subsequently the 
likelihood for considering the Q2 direction as a candidate for 
the next step is enhanced by a factor of nine (solid curve) 
which is seen to compensate exactly for the reduced lattice 
spacing. Also shown are the experimental data [251 . l2q | . 



In the Metropolis procedure )24| . a proposed step is 
always accepted if the associated energy change AU is 
negative, whereas it is accepted only with the probability 
exp(— AU/T) otherwise. The probabilities for accepting 
the reverse step are then exp(AC//T) or unity, respec- 
tively. So the ratio between the forward and backward ac- 
ceptance probabilities is P+/P- = cxp(— AU/T). Thus a 
continuous random walk of the type (f20|) can be treated 
approximately by means of a discrete random walk based 
on the Metropolis procedure. 

The treatment reported in Ref. 1] employed such a dis- 
crete random walk on the 5D lattice of potential energies. 
A bias potential of the form Vt>i as (Q2) = V0Q2/QQ was 
added to disfavor compact shapes and thus help to guide 
the walk towards the scission region, thereby speeding up 
the calculation; the strength used was checked to be suf- 
ficiently small that no further reduction would affect the 
mass yields. Thus the only physical parameter was the 
critical neck radius, cq, at which it was assumed that no 
further change in mass asymmetry would occur; typically 
the calculated mass yields were relatively insensitive to 
cq in the range of 2-3 fm and cq = 2.5 fm was adopted as 



the standard value in Ref. [I| ■ 

We are now in a position to ascertain the inaccuracies 
that are inherent in such an approach due to the specific 
lattice employed in Ref. [2] . We first examine the effect of 
the finite magnitude of the lattice spacing by comparing 
the results of the Metropolis walk with those of the cor- 
responding continuous process. This is illustrated in Fig. 
[1] for two typical cases (the others considered show com- 
parable effects). The continuous process and its discrete 
approximation tend to yield rather similar results for the 
fragment mass distribution. But there are some notice- 
able differences in the region of the asymmetric peak for 
the thorium isotopes considered. Such discrepancies sug- 
gest that the results are quite sensitive to the underlying 
potential energy surface in that particular region of the 
shape space (c/. the effect of modifying the Wigner term 
discussed in Sect. HVCj) . a feature that might help to 
achieve a better understanding of the potential energy. 

Secondly, we address the important role played by the 
relative size of the site spacings in the different lattice 
directions. The Metropolis walks carried out in Ref. 
1] treated all the five lattice directions equally which 
amounts to implicitly assuming that the underlying mo- 
bility tensor is isotropic, i.e. the mobilities are the same 
in all the lattice directions, fj, ~ I. Each lattice direc- 
tion has then an even chance for being considered as a 
candidate for the next Metropolis step. However, as the 
above analysis brings out, a change in the lattice spacing 
in one direction modifies the effective mobility in that 
particular direction: with an increased density of lattice 
sites, it takes a larger number of elementary Metropolis 
steps to go a given distance. Therefore, if the density 
of lattice sites is increased by a factor of N in the di- 
rection i then the corresponding mobility coefficient is 
decreased by the factor N 2 , pa — > m/N 2 . Conversely, 
to ensure that the evolution of the transport process re- 
mains unaffected by the increased density of lattice sites, 
the likelihood that the affected direction is being consid- 
ered as a candidate for the next step must be increased 
by that factor. For example, if the density of lattice sites 
in a particular direction is doubled, N = 2, then the like- 
lihood for considering that direction should be increased 
from i to 2 2 /(4 + 4) = \ to achieve the same evolution. 

These features are illustrated in Fig. [2] for the same 
two cases that were shown in Fig. Q] We consider the 
effect of changing the lattice spacing in the / direction 
which corresponds to the overall elongation as quantified 
by the quadrupole moment Q2 ■ When the lattice spacing 
in the / direction is decreased by a factor of three and 
the Metropolis walk is repeated with no other change, 
then the resulting fragment mass distribution is affected 
quantitatively, becoming noticeably narrower. However, 
as suggested by the analysis above, when the correspond- 
ing mobility is also increased by a factor of nine (by favor- 
ing the consideration of the / direction correspondingly) 
then the mass distribution reverts to its original form. 

The relative lattice spacings are thus intimately re- 
lated to the anisotropy of the effective mobility tensor. 



This basic feature makes the Metropolis walks performed 
in Ref. [l| seem somewhat arbitrary, because a differ- 
ent choice of lattice spacing, without any compensating 
change in the mobility coefficient, would generally lead 
to a different final result. It is therefore important to 
employ a mobility tensor based on a physically plausible 
form of the dissipation and to investigate the sensitiv- 
ity of the calculated mass distributions to that specific 
structure. We now turn to this central issue. 



III. INCLUSION OF DISSIPATION 

As discussed above, the simple random walk intro- 
duced in Ref. MJ] is most easily justified if the dissipa- 
tion tensor is isotropic in the employed lattice variables. 
Since such an idealized scenario is not likely to be realis- 
tic, we wish to study the effect of using a more plausible 
dissipation tensor, which is generally anisotropic and has 
a structure that varies from one shape to another. 

The potential energy of deformation, U(x), was calcu- 
lated [2j on a lattice of shapes introduced by Nix [27], [28j . 
As described in Appendix [B] each shape is composed 
of three smoothly joined quadratic surfaces. These 3QS 
shapes are characterized by the parameters q = {q^}- 
While these are in principle known functions of the lat- 
tice shape variables x — {Xi}> they are readily available 
only at the discrete lattice sites x — {I, J,K 7 L, M). 

The dissipation tensor can be determined from the 
rate of energy dissipation Q, which is a positive definite 
quadratic form in the shape velocities, 



Q = ^XiliAx)Xj = ^i^M)iu 



(24) 



The first term expresses Q in terms of the lattice vari- 
ables {xi}, while the second term uses the 3QS parame- 
ters {g M } so g^ v {q) is the friction tensor with respect to 
these variables (see Appendix [Bjl . Once g^iq) has been 
calculated (see below) , we may obtain the required dissi- 
pation tensor 'y(x) by the appropriate transformation, 






(25) 



We wish to determine -y(x) a t the various lattice sites, 
at which the parameters {x%} have integer values and we 
approximate the derivatives dq v /dxn in terms of differ- 
ences between the values of q v at the neighboring sites. 
Although this is a relatively rough approximation be- 
cause the dependence of q on x is generally not linear, it 
will sufifice for our present explorative purposes. 

In order to calculate g^iq), we thus need to know the 
dissipation rate Q(q,q). For that we employ here the 
"wall formula" for the one-body dissipation mechanism 
[7j. The underlying mechanism is the reflection of indi- 
vidual nucleons off the moving surface which generates a 
dissipative force that is rather strong due to the nucleonic 
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FIG. 3: The charge yields for neutron-induced fission of Pu 
and 236 - 234 U calculated with increasingly isotropic mobility 
tensors as obtained by using Eq. (|27p with / = 0.2, 1, oo (the 
latter is fully isotropic and has been indicated as / ;§> I), 
together with the experimental data [25l . l33| . Those in the top 
three panels are for (n t h,f) reactions j25j], while the data in the 
bottom panel is for (7,f) reactions leading to _B*« 8— 14 MeV; 
they include contamination from multi-chance fission [26J ]. 



Fermi motion. Because the individual nucleons reach the 
moving surface at random times, and from random di- 
rections, the associated force on the surface is stochastic, 
in accordance with the fluctuation-dissipation theorem 
[23j . While this idealization may not give a quantitative 
account of the actual dissipation rate in a real fissioning 
nucleus, it does serve well as a means to provide us with 
a mobility tensor that has a quasi-realistic structure. 

In its simplest form, the one-body dissipation rate in 
a deforming nucleus is given by the simple wall formula, 
Q = mpov § n(a) 2 d 2 a, where to, po, and v are the nu- 
cleonic mass, density, and mean speed in the interior, 
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FIG. 4: The charge yield for thorium isotopes calculated with 
increasingly isotropic mobility tensors as obtained by using 
Eq. (|27p with / = 0.2, l,oo (the latter being fully isotropic), 
together with the experimental data obtained with (7,f) reac- 
tions leading to E*~ 8 — 14 MeV; they include contamination 
from multi-chance fission 261. 



while h(a) is the normal surface velocity at the location 
a 0, [29| . It is elementary to show Q that the elements 
of the associated dissipation tensor are given by 



rPo« 



dp 2 (z) dp 2 (z) 
dqfj. dq v 



oz 



dz, (26) 



where p(z; q) is the transverse extension of the nucleus at 
the position z along the symmetry axis, for the specified 
values of the 3QS parameters q. The required quantities, 
namely p 2 (z) and its derivatives with respect to both z 
and the shape variables {g„}, can be expressed analyti- 
cally for the 3QS shape family J27| and so it is possible 
to calculate the elements of g for each specified shape. 



However, as explained in Appendix[B] certain elements 
of the dissipation tensor g may occasionally tend to zero 
(this happens when one of the three quadratic surfaces 
covers a negligible z interval so that this shrinking section 
ceases to contribute to the dissipation). Of course, the 
corresponding derivatives dq^/dxi diverge at the same 
time so the resulting elements of 7 remain well behaved. 
But, because it is impractical, for the time being, to 
calculate those derivatives with high accuracy, the cal- 
culation of 7 is correspondingly inaccurate, with some 
eigenvalues occasionally becoming unrealistically small. 

Fortunately, our main purpose here is merely to study 
the sensitivity of the calculated fragment mass distribu- 
tions to the structure of the dissipation tensor and we 
therefore perform the following isotropization procedure. 
If 7 is the original tensor, calculated as described above 
but renormalized so that its five eigenvalues {7™} are 
one on average, i.e. J2 n 7« = 5, then we define a more 
isotropic tensor 7 by modifying the eigenvalues, 



in 



[in + /]/(! + /) , 



(27) 



where the isotropization coefficient / is a positive num- 
ber. We see that the original friction tensor is recovered 
when / tends to zero, while it approaches isotropy when 
/ grows large. The corresponding modified mobility ten- 
sor p.(f) is then the inverse of'y(f). 

The sensitivity of the calculated charge yields to the 
degree of structure in the mobility tensor is illustrated 
in Figs. [SHU for the cases presented in Ref. [l|. In ad- 
dition to the experimental data, which are shown for 
reference, each plot shows the result of three different 
mobility scenarios: the idealized scenario (labeled f^>l) 
where the mobility tensor is isotropic, an intermediate 
scenario (/ = 1) in which the dissipation tensor is the 
average of the one calculated with the wall formula as 
described above and the corresponding directional aver- 
age, and a more structured scenario (/ = 0.2) in which 
the isotropic admixture is only 20 %. 

For the three neutron-induced cases, 239 Pu(n,f) and 
235 ' 233 U(n,f), the change in Y(Z{) is very small as one 
form of the mobility tensor is replaced by another, the 
most noticeable change being a slight narrowing of the 
asymmetric peaks for 233 U(n,f). The photon-induced re- 
actions, which are all calculated for E*=ll MeV, display 
a somewhat larger sensitivity. Generally, as the idealized 
isotropic mobility tensor grows more anisotropic there is 
a tendency for the symmetric yield component to become 
more prominent, but the quantitative effect is relatively 
modest. It is particularly noteworthy that the evolution 
from a symmetric yield for 222 Th to a mixed but pre- 
dominantly asymmetric yield for 228 Th remains present 
in all the scenarios. When comparing with these data, it 
should be kept in mind that they represent a range of ex- 
citations (E* f=a 8 — 14 MeV) and also are contaminated 
by second- and third-chance fission. 



IV. DISCUSSION 



B. Saddle shape 



We now discuss a number of interesting aspects that 
can be elucidated with the present treatment. 

The choice of shape degrees of freedom made in Ref. 
[2| , and of the specific 5D shape lattice used for the cal- 
culations of the potential energy, was guided in large 
part by physical intuition (using the somewhat vague 
but reasonable criterion that the typical energy change 
between neighboring sites should be of comparable mag- 
nitude). Our present studies suggest that this lattice 
of nuclear shapes was indeed well chosen, because the 
simple Metropolis walks provide mass yields that arc 
changed only moderately when a more refined treatment 
is made. Because the actual mobility tensor is not yet 
well known, it would seem prudent to employ a range of 
mobility scenarios. The spread among the results might 
then be taken as a rough indication of the uncertainty in 
the prediction. 

Furthermore, on the basis of our studies, it appears 
that the Metropolis walk, which is significantly faster 
than the Smoluchowski simulation (by 1-2 orders of mag- 
nitude), offers a very quick and easy means for obtaining 
practically useful fission-fragment mass distributions. 



A. Shape family 

We start by illustrating the importance of employing a 
shape family that has a sufficient degree of flexibility. For 
that purpose, we construct three-dimensional potential- 
energy surfaces by minimizing the full five-dimensional 
3QS surfaces with respect to the deformations of the two 
spheroids, en and €f2 (corresponding to the lattice indices 
K and L). Thus the shapes in the lower-dimensional 
space are characterized by only their overall elongation 
(represented by the lattice index /), their constriction 
(represented by the lattice index J), and the degree of 
reflection asymmetry (represented by the lattice index 
M). 

Figures [5] and [5] show the resulting charge distribu- 
tions for the cases presented in Ref. [1], together with 
the experimental data and our standard results based on 
the full 5D 3QS shape family. We see that the although 
the 3D calculations occasionally reproduce the qualita- 
tive appearance of Y(Z{) reasonably well, the reproduc- 
tion of the experimental data is generally far inferior to 
the results obtained with the 5D shape family. 

These examples demonstrate that it is important to 
employ a sufficiently rich family of shapes when seek- 
ing to describe the shape evolution during nuclear fis- 
sion. To obtain a reasonably flexible family of shapes, at 
least five shape degrees of freedom appear to be required, 
namely overall elongation, constriction, reflection asym- 
metry, and deformations of the individual prefragments. 



Because our calculational method emulates the actual 
equilibration process it is possible to gain insight into the 
shape evolution during fission. A particularly instructive 
case is presented by 222 Th: Our calculations yield a sym- 
metric mass distribution, in agreement with the experi- 
mental data, even though they are based on a potential- 
energy surface whose fission saddle point corresponds to 
a nuclear shape that is reflection asymmetric. 

The most relevant shapes are shown in Fig. [7J As 
the shape evolves from that of the ground state, it tends 
to pass near by the second (isomeric) minimum and the 
nucleus will typically remain trapped in that minimum 
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FIG. 5: (Color online.) Charge yields resulting from Metropo- 
lis walks on the 5D potential-energy surfaces associated with 
the 3QS shape family, together with the corresponding re- 
sults obtained with 3D surfaces generated by minimizing the 
5D surfaces w.r.t. the individual fragment deformations, en 
and 6f2. Also shown are the experimental data [251. |29] . 




Second minimum 



30 



40 50 60 

Fragment charge number Z f 




Outer saddle 




Most probable fragment mass division 
1 Re 





FIG. 7: (Color online.) Three shapes relevant for fission of 
222 Th: the outer potential-energy minimum (which is reflec- 
tion symmetric) where the fission isomeric state resides (top), 
the outer saddle which is asymmetric (middle), and the most 
probable fragment division which is symmetric (bottom). 



traverse the barrier with an asymmetry close to that of 
the saddle shape, P(Zf, /saddle) is concentrated around 
that asymmetry. However, beyond the saddle the pre- 
ferred asymmetry tends to become smaller (/ = 24) and 
the asymmetry eventually becomes peaked at symmetry 
(/ = 28). As the shape evolves further towards scission 
(/ = 32), P(Z{\ I) developes a minor asymmetric compo- 



FIG. 6: (Color online.) Similar to Fig. [5l but for fission of 
the isotopes 222,224,226,228 Th; the data are frQm Rgf j^ 



for quite some time before escaping, either back to the 
ground-state region or towards scission. (For that reason, 
we usually start our calculations at the second minimum, 
which reduces the required computational effort very sig- 
nificantly; we have of course checked that this does not 
alter the results.) As the figure shows, the outer min- 
imum of 222 Th is reflection symmetric while the outer 
saddle lies in a region of significant asymmetry Never- 
theless, the shapes evolve in such a manner that the final 
fragment mass distribution is centered around symmetry. 
More detailed insight into this evolution can be ob- 
tained by considering the charge-asymmetry distribution 
at specified values of the lattice index / which is a mea- 
sure of the overall quadrupole moment of the fissioning 
shape, P(Z{;I). This conditional distribution is shown 
in Fig. [5] for increasing elongations, starting at the value 
associated with the second saddle point, /saddle = 20. 
Because it is energetically favorable for the system to 
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FIG. 8: (Color online.) The evolution of the mass-asymmetry 
distribution from the region of the second saddle (/ = 20) 
towards scission is illustrated by the charge distribution 
P(Zf ; /) for increasing values of / (the lattice index giving the 
quadrupole moment Q2 of the nuclear shape) , as obtained for 
an ensemble of 10,000 Metropolis walks. 
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nent that presumably reflects the detailed (and possibly 
inaccurate) structure of the potential-energy surface in 
the scission region (see the discussion of the Wigner term 
in Sect. [JVC]). 

This result invalidates the commonly made assumption 
(see e.g. Refs. [30H33J ) that the character of the mass 
distribution, whether symmetric or asymmetric, is deter- 
mined by the character of the saddle shape. In contrast, 
analyzes of the type illustrated in Fig. [5] suggest that 
the structure of the potential-energy landscape in the 
entire region between the isomeric minimum and scission 
plays a role in determining the fragment mass distribu- 
tion. Obviously, any plausible model of the mass yields 
must take this into account. 



C. Wigner term 

The generality of our treatment makes it possible to 
exploit the remaining differences between calculated re- 
sults and experimental data to gain novel insight into as- 
pects of the fission process that would not otherwise be 
readily accessible. As an example, we consider here the 
shape dependence of the Wigner term in the macroscopic 
nuclear energy [34I439J . 

As mentioned above, the potential energies of Ref. [2(| 
were calculated as the sum of a macroscopic term and 
a microscopic (shell) correction, both being shape de- 
pendent. The macroscopic nuclear energy contains the 
so-called Wigner term proportional to \N — Z\. Its pres- 
ence is clearly visible in the systematics of nuclear masses 
which exhibit a V shape near isosymmetry, but its mi- 
croscopic origin is still not well understood, so it is nor- 
mally modeled as a phenomenological macroscopic term. 
In commonly employed models for nuclear masses [40l.l4l| 
the Wigner term is usually introduced without a shape 
dependence. However, such an ansatz presents a signifi- 
cant (but often ignored) problem in fission where a sin- 
gle nucleus is transformed into two separate nuclei, each 
having its own Wigner term. (That each fragment nu- 
cleus must give a Wigner contribution similar to that of 
the original nucleus is evident from the phenomenolog- 
ical form of the term [3a, |2Z|.) Thus the Wigner term 
must double in magnitude during the evolution from a 
single compound system to two separated fragments, so, 
consequently, it must depend on the nuclear shape. 

The calculations of fission barriers in for example Refs. 
@, IH, |36| have since 1989 employed a postulated shape 
dependence that relates the increase of the Wigner term 
to the decrease in the amount of communication be- 
tween the two fragments due to the shrinking of the neck. 
The resulting Wigner term changes gradually as the nu- 
clear shape evolves and affects the potential-energy land- 
scape correspondingly (an illustrative figure was given in 
Ref. |3J|)- However, it may be argued that the change in 
the character of the fissioning system from mononuclear 
to dinuclear occurs more abruptly than implied by the 
currently prescribed shape dependence. To elucidate the 
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FIG. 9: (Color online.) The sensitivity of the charge yields 
to the shape dependence of the Wigner term is illustrated 
for 226 Th, using either the standard potential-energy surface 
2] , in which the Wigner term changes gradually as the shape 
evolves, or a modified energy surface obtained with a Wigner 
term that remains constant until scission occurs; the three 
panels show results calculated with increasingly isotropic mo- 
bility tensors as obtained by using Eq. (|27[) with / — 0.2, 1, oo; 
the experimental data are also shown [261 ]. 



importance of the shape dependence of the Wigner term, 
we have considered an alternate form in which the term 
is kept constant up to the scission configuration, i.e. un- 
til the neck radius has shrunk below the specified critical 
value Co- 

In Fig. [9] we illustrate how such a modification of the 
calculated potential-energy surface affects the calculated 
charge distribution of 226 Th, for which the impact is par- 
ticularly noticeable. There are two significant differences 
between the results of the two sets of calculations. One 
is a change in the relative importance of symmetric and 
asymmetric fission, with the constant Wigner term lead- 
ing to more asymmetric yield. The other is a shift in the 
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FIG. 10: The charge yield for the 234 U(7,f) reaction extracted 
from Metropolis walks using either the macroscopic expres- 
sion for the level density, a a = A/eo, with the standard value 
eo=8 MeV and also eo=10 MeV, or the microscopic expression 
(|28|) ; the experimental data are included [2a ]. 



location of the asymmetric yield peaks, from being sev- 
eral units on the outside of the observed values towards 
a better agreement with the data. Both effects depend 
significantly on the structure of the mobility tensor and 
could, in principle, be of help in discriminating between 
different models of the dissipation. 
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pression ([28]) . On the whole, the calculated mass yields 
are remarkably unaffected by this change, probably be- 
cause each random walk visits quite a large number of 
lattice sites, so the shell effect, which tends to fluctuate, 
largely averages out. The most noticeable effect occurs 
for the 234 U(7,f) reaction where the use of the micro- 
scopic expresison leads to more yield in the symmetric 
region (where our standard calculation underpredicts the 
yield). This case is shown in Fig. [TOl where we have also 
included the result of using the macroscopic formula with 
a larger value of eo- Due to the relation E* = ciaT 2 , this 
change causes the local temperature to increase and, as a 
result, the symmetric yield is increased. As it turns out 
the effect of making this increase in eo in the macroscopic 
formula is practically identical to the effect of replacing 

^raacro oy Ctniicro- 

We wish to point out that our present studies do not 
consider the effect of pairing which may be included in a 
simple approximate manner by backshifting the excita- 
tion energy E*(x) by the pairing gap A(x) when calcu- 
lating the local temperature T(x)- Such an undertaking 
would require knowledge of the shape-dependent pairing 
gap, A(x). While this quantity was of course calculated 
at each lattice site when the potential-energy surfaces 
were generated Q, it was not tabulated seperately, so it 
is not presently available. We must therefore leave this 
interesting issue for future study. 



E. Scission model 



D. Level density 

The microscopic part of the potential energy, SU s h (x) , 
is due to the deformation-dependent variations in the 
single-particle level densities in the effective field of the 
fissioning nucleus. This structure also affects the depen- 
dence of the nuclear temperature T on the excitation 
energy E* . Because a change of the local temperature 
affects the local diffusion rate but not the drift rate, a 
change in T(E*;x) mav influence the evolution of the 
shape distribution P{x)- 

In order to explore the importance of this effect, we 
replace the standard Fermi-gas level-density parameter 
cia = A/eo m the formula E* = ciaT 2 by a "shell- 
corrected" generalization suggested by Ignatyuk [42j , 



a A (E*;x) = o-A 



l + (l-c- E */ £ — ) WshW 



E* 



, (28) 



where -Edamp characterizes the gradual dissolution of the 
shell effects as the excitation energy is increased. Wc 
shall use e = 8 MeV and E damp = 18.5 MeV. At low 
excitation, E* — > 0, a a tends to a[l + <5£/ s h/-Edamp] ) while 
it approaches cla monotonically as E* is increased (when 
E* ^> -Edamp the exponential is close to zero and also 
SU sh /E* < E d&mp /E* < 1 since \SU sh \ < E damp ). 

We have examined the effect of replacing the macro- 
scopic formula cia — A/eo by the above microscopic ex- 



Finally, we wish to illustrate the importance of the pre- 
scission shape evolution by comparing with the statistical 
scission model [4j, [44| . It assumes that the scission con- 
figurations, which are here those having c(x) — c 0j are 
populated in proportion to exp(— U(x)/T(x))- 

Such calculations are, by definition, entirely static in 
nature and do not include any dynamical effects. While 
scission models often yield reasonable agreement with the 
observed mass distributions, they are not universally suc- 
cessful and their failures suggests that mass splits are sen- 
sitive to the pre-scission evolution. In particular, some 
scission configurations may not be easily reached due 
to the presence of intermediate ridges in the potential- 
energy landscape. Such is presumably the case, for ex- 
ample, in the recently reported fission of 180 Hg in which 
symmetric splits are strongly favored by the energetics 
in the exit channel but are most likely prevented dynam- 
ically by the presence of a potential-energy ridge, thus 
causing the fragment mass distribution to become asym- 
metric mm- 

But even in cases where there are no such prominent 
obstacles present in the potential-energy landscape, the 
statistical model is often inaccurate. This is illustrated in 
Fig.[TT]in which its results for the plutonium and uranium 
cases are compared with both the experimental data and 
the results obtained with the simple Metropolis walk in- 
troduced in Ref. |1|. The statistical scission model re- 
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FIG. 11: Results of the scission model (see text) for the same 
plutonium and uranium isotopes as shown in Fig. 6, together 
with the corresponding experimental data and the results of 
the simple Metropolis walk introduced in Ref. [j|. 



produces the qualitative appearance of Y(Z{) reasonably 
well, but the quality of the correspondance with the data 
is inferior to that provided by the Metropolis transport 
calculation. 



V. CONCLUDING REMARKS 

Generally, calculations of the type discussed here per- 
tain to the idealized limit of stongly damped motion 
where the inertia plays no role. The evolution of the nu- 
clear shape is then akin to Brownian motion. Whether 
in fact the shape evolution during fission has such a char- 
acter is still an open question, but we believe that this 
simple limit provides a useful reference scenario. 

In a first exploration of this physical picture, it 



was recently shown that simple random walks on five- 
dimensional potential-energy surfaces lead to remark- 
ably good agreement with experimental data on fission- 
fragment mass distributions |l|. We have examined this 
method in depth, studying the importance of a number 
of effects that could be expected to influence the results. 

First we generalized the discrete random walk on the 
fixed lattice sites where the potential energy is available 
to a continuous diffusion process, thus enabling an as- 
sessment of the importance of the finite lattice spacing. 
We found that there is typically little difference between 
the results of the two treatments, although we observed 
some deviations in certain limited regions. 

Such a simple random walk, whether discrete or con- 
tinuous, is physically reasonable only if the dissipation 
tensor happens to be isotropic in the particular shape 
variables employed, which is generally not expected to 
be the case. This important feature was illustrated by 
studying the effect of inserting additional sites between 
the standard lattice sites which reduces the mobility in 
the affected direction. 

The main objective of the present study is therefore 
to elucidate how the results of the idealized treatment 
may change when a more realistic dissipation tensor is 
employed. For this purpose we employed the dissipa- 
tion tensor suggested by the simple one-body dissipation 
wall formula and introduced an isotropization procedure 
which allowed us to examine a continuum of scenarios 
ranging from perfect isotropy (corresponding to the ide- 
alized random walks discussed above) to the relatively 
large degree of anisotropy displayed by the calculated 
friction tensors. Generally we found that the resulting 
fragment mass distributions are rather insensitive to the 
degree of anisotropy, except, in certain cases only, for 
rather extreme anisotropies that are probably unrealistic 
and may arise from certain numerical problems. 

We then examined a number of additional relevant as- 
pects. First we demonstrated the importance of using 
a sufficiently rich family of shapes by comparing results 
based on the full 5D potential energy surfaces with anal- 
ogous calculations in a reduced 3D deformation space ob- 
tained by constrained minimization. The importance of 
using sufficiently flexible shapes was borne out in partic- 
ular by the isotope 222 Th whose mass yield is symmetric 
even though the outer saddle shape is asymmetric. 

As an example of the potential of the present method 
for elucidating novel aspects of nuclear properties, we 
examined the sensitivity of the calculated mass distribu- 
tion to the character of the deformation dependence of 
the Wigner term in the macroscopic energy functional. 
Specifically, it was demonstrated for 226 Th that when 
the previously employed rather gradual shape depen- 
dence was replaced by a more abrupt transition from the 
mononucleus to the dinucleus then the symmetric compo- 
nent of Y(Z{) is reduced and its asymmetric peaks move 
towards better agreement with the experimental data. 

We also examined the effect of taking account of the 
deformation-dependent bunchings in the single-particle 
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level densities which affect the energy dependence of the 
temperature for a given shape and thus the local diffu- 
sion rate relative to the drift. The calculated mass yields 
are relatively unaffected by this refinement, except for a 
favorable shift in the amount of nearly symmetric mass 
splits, but it should be kept in mind that the effect of 
the (presently unavailable) shape-dependent paring gap 
is still left for future study. 

Finally, comparisons with mass distributions arising 
from a statistical population of the scission configurations 
demonstrated the importance of pre-fission shape evolu- 
tion. Typically, our transport calculations yielded signif- 
icantly better agreement with the experimental data. 

For some of the cases considered here, most notably 
239 Pu(n, f), the calculations reproduce the experimen- 
tally observed yields very well and are rather robust with 
respect to model variations, such as changes in the mo- 
bility tensor. Hence, for the comparisions with data to 
be informative, it appears to be important to include also 
cases that exhibit larger sensitivity, such as the thorium 
isotopes. 

In conclusion, our studies suggest that the sim- 
ple Metropolis walk [1] on the previously calculated 
potential-energy lattice |2j indeed presents a useful 
calculational tool for obtaining the approximate form of 
fission-fragment mass distributions for a large range of 
nuclei. For more accurate results it is necessary to invoke 
also the dissipative features of the shape evolution as 
represented by the shape-dependent mobility tensor. 
The shape evolution then resembles Brownian motion 
in an anisotropic (and non-uniform) medium. However, 
because the dissipation mechanism is not yet as well 
understood as the potential energy, we propose to 
make a series of calculations with mobility tensors that 
display different degrees of anisotropy and then use the 
ensuing spread in the results as an indication of the 
uncertainty of the predicted mass yield. The results 
obtained in this manner are often remarkably robust. 
Consequently, the method may be of practical use for 
calculating fission-fragment mass distributions for any 
of the thousands of nuclei for which the required 5D 
potential-energy surface is already available. 
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Appendix A: Lattice interpolation 

The nuclear shapes are characterized by the five shape 
parameters {x n } = (%i, X2, X3, X4,Xs) which are collec- 
tively denoted by X- But the potential energy of deforma- 
tion, C/(x), is known only on a five-dimensional Cartesian 
lattice on which the shape parameters take on the inte- 
ger values X = {X n } = (/, J, K, L, M), whose ranges are 
given in Ref. Q. (We use X to denote a capital X-) We 
describe here how the potential energy for arbitrary x 
values can be obtained by pentalinear interpolation, i.e. 
an interpolation scheme that yields a function U({x n }) 
that is linear in each of the five variables Xn inside each el- 
ementary hypercube. [The resulting function is identical 
to the Taylor expansion around the "lower-left" hyper- 
cube corner, X = (7~, J~ ,K~ ,L~ ,M~), keeping only 
terms of first order in each variable and approximating 
all derivatives by the corresponding central differences.] 

We assume that the given shape parameter x h es 
within the domain covered by the lattice and start by 
identifying the surrounding elementary hypercube. Its 32 
corners are given by the indices (I*, J ± , A ± , L ± ,M ± ), 
where X~ = [x„] (i.e. the integer part of Xn) and 
X+ = X- + 1 so X- < Xn < X+ for n = 1, . . . , 5. 
In a single dimension, the interpolated value would be 

U(x) = U(X-)(X+- X ) + U(X+)( X -X-), (Al) 
so we may readily generalize to five dimensions, 

U(X) = ]T ijklm 17(7"* J- J , K-\ L-\ A<r m ) (A2) 

ijklm 

x {P- Xi)(J J -X2)(K k - X3)(L l - Xi )(M m - xn), 

where the summation indices each take on the values =f1. 
It is easy to verify that this is indeed correct: Within the 
local hypercube (within which x i s located) the above 
expression is linear in each of the five x variables and 
it yields the correct matching because when x coincides 
with a lattice site, x = X = (/, J, K, L, M), then I~ = I, 
..., M ~ = M so X~ — Xn = and X+ — Xn = 1, so only 
the term with i = +1,. .. ,m = +1 contributes, yielding 
U(x) = U(I,J,K,L,M). 

The driving force F = — Vt/(x) can be obtained by 
taking the derivative of the above expression (|A2|) . Thus 

Fi(x) = -Yl V klm U ^ J ~ J < K ~ k > L ~ l > M ~ m ) 

ijklm 

x (J*- X 2)(K k - X3 )(L l - Xi )(M m - X5 ), (A3) 

and analogously for the other four directions. Thus, 
within the local hypercube, the force component F n does 
not depend on x„, as is consistent with the fact that the 
potential is locally linear in %„. 

A similar scheme is used to calculate other quantities 
for arbitrary shapes, such as the neck radius c(x) an d 
the dissipation tensor ~f(x)- 
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Appendix B: 3QS shape family 

The three-quadratic-surface shape family introduced 
by Nix [28| consists of axially symmetric shapes for which 
the square of the local radial distance p(z) to the surface 
is given by three smoothly joined quadratic surfaces, 

( a\ - (a 2 J4){z - h) 2 , h~c x <z< Zl , 
p 2 (z) = 1 a 2 3 - (a 2 3 /4)(z - £ 3 ) 2 , z x <z< z 2 , (Bl) 
{ a 2 - (a 2 2 /c 2 )(z - £ 2 ) 2 , z 2 <z<£ 2 + c 2 . 

Thus nine numbers are required to specify the nuclear 
surface. Two of these are eliminated due to the continuity 
of p{z) and its derivative at z\ and z 2 , and the length 
parameter u = [\{a\ + a?,)] 1 / 2 governs the overall scale. 
The remaining six numbers are then determined by six 
dimensionless shape parameters {q u }, 



o\ = 






Oil 



h 



a 2 = 



a 3 = 



(B3) 



Furthermore, if the shapes are required to have a given 
center of mass, then the parameter ai is determined once 
the other five have been specified. All the shapes in the 
IJKLM lattice have their center of mass at the ori- 
gin, which effectively reduces the six-dimensional {q v } 
space to the five-dimensional shape space covered by the 
IJKLM lattice. 

For the evaluation of the dissipation tensor 7^ we need 
the derivatives dp 2 /dq v as well as dp 2 jdz which, though 
somewhat involved, can be expressed analytically [271 ] . 

The 3QS family includes shapes for which one of the 
three segments covers only a negligible z interval. The 
contribution from such a segment to the dissipation rate 
is then also negligible and, as a result, so are the associ- 
ated elements of the dissipation tensor. These singulari- 
ties are numerically inconvenient and must be addressed. 
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